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Abstract 



o , 

, We discuss the relationship between kinetic equations of the Fokker- 

Planck type (two linear and one non-linear) and the Kolmogorov (a.k.a. 
master) equations of certain A^-body diffusion processes, in the context 
of Kac's propagation-of-chaos limit. The linear Fokker-Planck equa- 
tions are well-known, but here they are derived as a limit ^ oo of 
' a simple linear diffusion equation on 3A^ — C-dimensional A^-velocity 

spheres of radius oc ^/N (where C = 1 or 4 depending on whether 
the system conserves energy only or energy and momentum). In this 
case, a spectral gap separating the zero eigenvalue from the positive 
spectrum of the Laplacian remains as ^ cxd, so that the exponen- 
5^ \ tial approach to equilibrium of the master evolution is passed on to the 

limiting Fokker-Planck evolution in M^. The non- linear Fokker-Planck 
equation is known as Landau's equation in the plasma physics literature. 
Its A^-particle master equation, originally introduced (in the 1950s) by 
Balescu and Prigogine (BP), is studied here on the 3A^ — 4-dimensional 
A^-velocity sphere. It is shown that the BP master equation represents 
a superposition of diffusion processes on certain two-dimensional sub- 
manifolds of determined by the conservation laws for two-particle 
collisions. The initial value problem for the BP master equation is 
proved to be well-posed and its solutions are shown to decay exponen- 
tially fast to equilibrium. However, the first non-zero eigenvalue of the 
BP operator is shown to vanish in the limit — > oo. This indicates that 
the exponentially fast approach to equilibrium may not be passed from 
the finite- A^ master equation on to Landau's nonlinear kinetic equation. 

(c)2004 The authors. This paper may be reproduced for noncommercial purposes. 
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1 Introduction 



Kinetic equations play a crucial role in the transport theory of gases and plas- 
mas, in particular for studying the approach to equilibrium. Apart from rigor- 
ous mathematical studies of their solvability properties and the classification 
and description of their solutions, it is essential to establish their validity. The 
validation of a kinetic equation consists in its derivation from some deeper, 
deterministic microscopic model, for instance from the classical Newtonian 
dynamics of an isolated system of many interacting (point) particles. Clearly, 
a complete validation automatically involves existence and uniqueness results 
for the kinetic equation that is being validated, and it also involves proving 
some version of the second law of thermodynamics. A priori knowledge of 
existence and uniqueness of solutions to the kinetic equation can aid the proof 
of its validity, while in the absence of such a priori knowledge the successful 
validation would yield existence and uniqueness for the kinetic equation as 
a corollary. Unfortunately, validation has turned out to be a very difficult 
problem. Even the probably most re-known and most studied of the kinetic 
equations, namely Boltzmann's equation, has been validated only in a few 
"simple" situations |Lan75[ |5^39T1 1( ^1 P94j . 

Meanwhile, inspired by the pioneering work of Kac |Kac56j a large body 
of literature has accumulated in which the deterministic N-body dynamics is 
replaced by an interacting stochastic Markov process which preserves, in each 
binary interaction, at least particle number and energy, but preferably also 
momentum and angular momentum, and which is designed to formally lead 
to the same kinetic equation that one expects from the "physical" N-body 
system in the infinitely many particles limit through a law of large numbers. 
Technically, the law of large numbers for the stochastic evolution of a family of 
individual systems of particles is equivalent to studying Kac's propagation 
of chaos limit N —>■ oo for the corresponding ensemble of individual systems, 
the evolution of which is being described by the Kolmogorov equation (called 
master equation in the physics literature) for the selected Markov process. Like 
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Liouville's equation, the Kolmogorov equation is a linear deterministic partial 
differential equation for the ensemble probability density on 6A^-dimensional 
phase space. Unlike Liouville's equation, the Kolmogorov equation typically 
defines a contraction semi-group instead of a group as does Liouville's equa- 
tion. Hence, the relaxation of the ensemble density to a uniform density is 
now built into the ensemble evolution, and because of Boltzmann's result that 
in the limit N oo almost every point of phase space corresponds to the 
Maxwellian velocity distribution, the second law of thermodynamics is a fore- 
gone conclusion. In this sense, the Kac approach is simpler than the original 
validation problem for kinetic equations, but it still retains some flavor of val- 
idation. From the perspective of validation, one could say that the Kac type 
approach goes "half the way" toward what one would like to prove. From 
the perspective of the mathematical analysis of the kinetic equations itself, 
the Kac type approach offers a new angle of attack to establish existence and 
uniqueness of the evolution and the relaxation to equilibrium in those cases 
where these results have not yet been obtained by direct PDE methods. 

Yet, the question whether the information obtained for a linear master 
(Kolmogorov) equation for finite N (such as global existence and uniqueness 
of solutions, as well as exponentially fast approach to equilibrium) carries on 
to the typically nonlinear kinetic equation which is expected to arise from the 
master equation in the limit — oo, has turned out to be more subtle than 
originally anticipated. The current state of the art of this approach for short 
range binary processes is presumably the work |CCL02j . 

Our primary concern in this paper is the master equation approach to 
certain kinetic equations that arise in the theory of systems with long range 
interactions, such as Coulomb plasmas and Newtonian gravitating systems. In 
particular, we discuss a master equation, originally introduced (in the 1950s) 
by Balescu and Prigogine (BP), that leads formally to the spatially homoge- 
neous Landau kinetic equation |Lan37j . which plays a fundamental role in the 
classical transport theory of Coulomb plasmas )Bal88| IHin83j . Considering 
here only the one-component case, the Landau equation for the particle den- 
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sity function f {.;{): M.'^ ^ M+ on velocity space at time t G M+ has the form 

dtf{v; t) = d^- [ Ql{v, w) ■ {d^ - d^) {f{v; t)f{w; t)) d'w, (1) 

where Ql{v, w) is the Landau colhsion kernel 

Qi^iv, w)=d^'\v-w\ = \v-w\-' Pt^ (2) 

with P^_w the projector onto the plane perpendicular to v — w E M^. We 
remark that with the help of the so-called Rosenbluth potentials |RMcD,T57] of 
/, the Landau equation can be recast as a nonlinear Fokker-Planck equation, 
the form that is better-known to the astrophysics community. Formally, the 
Landau equation satisfies the standard conservation laws of mass, momentum 
and energy; also the H-theorem holds. Thus, it is commonly believed that at 
late times the solution / evolves into the Maxwellian /m associated, via the 
conservation laws, with the initial data /q. Estimates of the relaxation time 
are usually obtained by linearization of the equations, but without estimates 
of the time it takes the dynamics to reach the linear regime. Despite its 
physical importance, the mathematically rigorous confirmation of the expected 
behavior of the solutions to this equation is lacking. Only very recently has the 
spatially homogeneous Landau equation attracted some attention in the PDE 
literature, where it has been studied as a member of a more general family of 
equations with kernels 

Q{v,w) = \v-w\'+^ Pt^ (7>-5) (3) 

7-5 

formally associated with 1/rT-i force laws. However, for Coulomb and Newton 
interactions (7 = —3), these PDE methods have provided only weak existence 
results |Vil96| IVil98aj . leaving the questions of uniqueness, regularity and ap- 
proach to equilibrium largely unanswered, except for initial conditions close 
to equilibrium |Guo02j or locally in time |Zha94j . 

As regards the Balescu-Prigogine A^-particle master equation for the Lan- 
dau equation, we will show that its initial value problem is well-posed, and 
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that its solutions approach equihbrium exponentially fast, on the 3N — 4- 
dimensional A^- velocity sphere of constant mass, energy, and momentum. How- 
ever, in the limit N ^ oo, with energy and momentum scaled so that the corre- 
sponding quantities per particle are constant (the A^-velocity sphere has radius 



This indicates that the exponentially fast approach to equilibrium described 
by the finite- A?" master equation may not be passed on to Landau's nonlinear 
kinetic equation. To resolve this issue requires a more detailed knowledge of 
the spectrum of the BP operator. 

While our efforts have not yet revealed all the details of the BP operator 
spectrum, we have discovered that the BP master equation represents a su- 
perposition of diffusion processes on certain two-dimensional sub-manifolds of 
IR^^ determined by the conservation laws for two-particle collisions. This has 
prompted us to study in more detail the completely solvable cases in which 
the Kolmogorov equation is just the linear diffusion equation dtF = AF on a 
3N — C-dimensional many-velocity sphere, where C = 1 or 4. The underlying 
stochastic processes are the perhaps simplest single + binary processes pre- 
serving either particle number N and total energy J^^^ | I'yfel^ = Ne {C — 1), 
or particle number, total energy and total momentum '^^^i — Nu (C = 4). 
The results are interesting enough to deserve being included in this paper. 

Thus, we show explicitly that in the limit N ^ oo we obtain an essentially 
linear Fokker-Planck equation for the particle density function f{-;t) : ^ 
M+ on velocity (=momentum) space at time t e M_|_, which for C = 4 has the 
form 



oc VN), the first non-zero eigenvalue of the BP operator is shown to vanish. 





for initial data f{v;0) > having mass per particle 




(5) 



momentum per particle 




(6) 
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and an energy per particle 

/ l-\v\'f{v;0)d'v = e, (7) 

which itself is a sum of the energy per particle in the center-of-mass frame, Eq, 
and the energy per particle of the center-of-mass motion, ecm = viz. e = 

Eq + It is easy to show that with such initial data the mass per particle 

f{v;t)d^v, the momentum per particle f^svf{v;t)d^v, and energy per 
particle J^^ ^Iv]"^ f{v;t)d^v are conserved during the evolution. Moreover, for 
such data the solution / of Q evolves, as t ^ oo, into the drifting Maxwellian 
/m associated, via the conservation laws, with the initial data f{v; 0); viz. 

and it does so exponentially fast and with monotonically increasing relative 
entropy 

S{f\fM) = - [ f{v;t)\n^p^d'v, (9) 

so that an H theorem holds. The treatment without momentum conservation 
(C = 1) is similar; of course, u does not show in this case. 

We remark that our derivation of (jH) together with ©, (0), and ((Tj) from 
an isolated system of particles preserving energy and momentum may not 
be new; however, since we could not find it in the literature, this interpretation 
of (jH) may perhaps not be so widely known. Indeed, (jH) is usually associated 
with an Ornstein-Uhlenbeck process for a swarm of individual, independent 
particles, the velocities (in M^) of which are being thermalized through contact 
with a heat bath of temperature T. In this case 2eq in is replaced by 3T 
and the restrictions (jH)) and ((Tj) on the initial data have to be dropped. 

In the remainder of the paper, we first set up the master equation approach 
for isolated spatially uniform systems in general. Then, to have a simple 
illustration of validation a la Kac, we first discuss the diffusion equations on 
the 3A^ — C-dimensional velocities spheres and derive the linear Fokker-Planck 
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equation(s) in the limit N ^ oo. Then we turn to the Balescu-Prigogine 
master equation and its putative N ^ oo hmit, the nonhnear Fokker-Planck 
equation due to Landau. 

2 Ensembles of isolated systems 
2.1 The velocity manifolds 

Let {Vq}^]^ denote an infinite ensemble of identically distributed random 
vectors taking values in M.^^ . Each vector V — {vi,...,Vn) G M.^^ represents 
a possible micro-state of an individual system of N particles with velocities 
f^i — {vii,Vi2,Vis) e R^. The positions of the particles are assumed to be 
uniformly distributed over either a periodic box or a container with reflecting 
boundaries and have been integrated out. The micro-state is assumed to evolve 
in time according to some stochastic process which conserves 

m{V) = N (mass of V) (10) 

and 

^('^) = ^ E l^^l' (energy of y), (11) 

where I^v = {1, N}; in the periodic box also 

p{V) = ^ Vk (momentum of V) (12) 

fceiiv 

is preserved. We are only interested in ensembles of particles systems in 
which all members have the same energy, or same energy and same momentum. 
Thus, depending on the circumstances, an ensemble consists of vectors Va 
taking values either in the 3N — 1-dimensional manifold of constant energy 

Mf -{V : ^ E l^^l' = ^ > O}' (13) 

or in the 3N — 4-dimensional manifold of constant energy and momentum 
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Each such manifold is invariant under the process which generates the evo- 
lution of an individual, isolated body system, which in turn traces out a 
trajectory on one of these manifolds M^^^-*^ or M^^~^. 

The manifold M^^~^ is identical to a 3A^ — 1-dimensional sphere S^f^^ 
of radius \/2Ne and centered at the origin of M?^] the manifold M^^~'^ is 



identical to a 3A^ — 4-dimensional sphere of radius y/2Neo and centered at 
U = {u, u), embedded in the 3(A^ — l)-dimensional affine linear subspace of 
^SAf gi^g^ t/+L3^-3, where L^^-^ = ^^^p { V G R^^ : J^keiN = O}. 
In the following, when we write , we mean S^^^^^ C L^^"^. 

2.2 Master equations 

Any ensemble {Vq,}^^ at time t is characterized by a probability density on 
either M^^^-*^ or M^^~^, for simplicity denoted F'^^\V] t). The time evolution 
of t) is determined by a master equation on n or 2? n 

2H 



If-'] 



where C'^^'^ is a positive semi-definite operator on £^(M^^~'*) or £^(M^^~^); 
(jl5p is the Kolmogorov equation adjoint to the underlying stochastic process. 
In general, the operator C^^^ has a non-degenerate smallest eigenvalue and 
corresponding eigenspace consisting of the constant functions. Since all par- 
ticles are of the same kind, we consider only operators C'^^'^ and probability 
densities F*^^^ which are invariant under the symmetric group applied to 
the N components in R? of V . Also, the density F'^^\V] t) has to satisfy the 
initial condition WmtiQ F^^\V]t) = fI^^\v). 

3 The diffusion master equation 

Since it is instructive to have some explicitly solvable examples, in this section 
we consider first the case of an isolated gas in a container, then we turn to the 
case of an isolated gas in a periodic box. 
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3.1 Gas in a container 
3.1.1 Finite N 



Taking 

£W = -A„3.-., (16) 

the master equation is then simply the diffusion equation on M^^^^ = S^j^^, 

a,FW(V;t) = A„3.-.FW(V;t). (17) 

Now, the spectrum and eigenfunctions of the Laplacian on a D-dimensional 
sphere are well-known. Since AgSiv-i = ^^Agsiv-i it will be enough to consider 
them on the unit sphere; the results can then be adapted to S^^^^ by simple 
scaling. With D = 3N — 1, the spectrum of — Agsiv-i is given by {Agg^jv-ij^lo 
with Ag3^jv-i = j(j + 3A^ — 2), and the eigenspace for Agi^jv.i is spanned by 
the restrictions to S'^^^^ C of the harmonic polynomials which are ho- 
mogeneous of degree j in M^^; when j > this restriction to S^^~^ C 
has to be non-constant. Solutions of the diffusion master equation which are 
invariant under the symmetry group Sn acting on V, however, can be ex- 
panded entirely in terms of eigenfunctions having that same symmetry. The 
simplest such eigenfunctions are the restriction to S^^'^ of the polynomials 
of the form Pj^\v) = "^^^i^ Pj{vk) where the pj{v)^s are harmonic polyno- 
mials which are homogeneous of degree j in (however, the special case of 
P2{v) = vf + V2 — 2f| simply leads to the constant function on S^^~^ and, 
hence, does not lead to an element of the eigenspace of Ag3iv-i)- In the next 
more complicated case the SAr-invariant eigenfunctions for j = ji + j2 are of 
the form P-^lj^iV) = Efceiiv ^ iei'^^l^Phi'"k)Pj2{vi), restricted to S^^-^; etc. 

ir3Af-l _ (s3Af-l 



Thus, the ^Ar-symmetric solutions of equation (|T7j) on Mf -1 = §^j^ are 
given by a generalized Fourier series 



where {G'y;^(V), £ G ©,} are the eigenfunctions of — AgSiv-i spanning the S^- 



symmetric eigen-subspace for ^^^sn^i for j G N, with 3j C N the set of indices 

\/2Ne 



/2Ne 



labeling the degeneracy of the j-th eigenvalue, and the Fj^^ are the expansion 
coefficients. (Although the eigenf unctions are quite explicitly computable, we 
here refrain from listing them; we shall only work with some simple eigen- 
functions for the purpose of illustration.) Evidently, the ensemble probability 
density function on S^j^^ evolves exponentially fast into the uniformly spread- 



out probability density 



which is the eigenfunction for A 



(0) 

3JV-1 
\/2Ne 



0. 



3.1.2 The limit ^ oo 



To discuss the limit N oo for the time-evolution of the ensemble, we con- 
sider the time-evolution of the hierarchy of n-velocity marginal distributions 
...,Vn]t) with domains {{vi, . . . , : ELi \^k\^ < 2^^^} C M^", 

which obtains by integrating (ITHl) over the available domains 

of the remaining N — n velocities variables, thus 



f n e 



i(i+3Af-2) 
2Ne 



t 



(19) 



where F^^}^\vi, . . . , i;„) and g^"'j^\vi, . . . , v„) are the corresponding integrals 



of 



V2Ne 



-1 



and G^^(^{V) over 



The n-velocity marginal of 

g3(7V-n.)-l 



1 

is 



3Af-l 
\/2Ne 



,Vn) 



33Af-l| 



2Ne 



2 



1 1 V^^*^ I I 



3(iV-n)-l 



(20) 



and converges in the limit — oo to the n-velocity Maxwellian 

n 

/r(«„...,«„) = n(4!i)^-p(-ii^^p) 



(21) 



k=l 



with domain M . Hence, one obtains the same equilibrium statistical mechan- 



ics from 



V2Ne 



viewed ClS db probability density on 



ir3Af-i 



as one obtains 



from the conventional Boltzmann-Gibbs micro-canonical equilibrium ensem- 
ble Z^^6{^ J2k&N ~ viewed as a stationary probability "density" on 
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^3N velocities of the perfect classical gas of particles at energy A^. 

Note that for finite A^, the results differ slightly. 

As for the partially integrated S'Tv-symmetric eigenfunctions of — AgSiv-i, 
suitably normalized they converge pointwise to a compatible family of Sn- 
symmetric functions on R'^", 

lim c/5"'^Hvi, ...,Vn)= gf}{vi, ...,Vn). (22) 

A detailed calculation, to be presented elsewhere, shows that the limit func- 
tions g^^J[v\^ . . . , i>„) are identically zero unless i belongs to a certain subset 
Dj C ©j. If £ G ©j, each g^^J{y\, . . . , i>„) turns out to be, for all n, one of the 
well-known jRisk96l| eigenfunctions of the Fokker-Planck equation in M^" (see 
(|26|) here below). Each eigenfunction is given by the n- velocity Maxwellian 
multiplied by a (symmetrized) product of Hermite polynomials, one in 
each component of the n velocities . . . , f„, of total degree j. 

Regarding the spectrum of — AgSiv-i, it is readily seen that in the limit 
N —>■ oo we have 



jini^liur„ = lir , (23) 



Thus, the whole spectrum remains discrete and, as is well known |Risk96j . 
coincides with the spectrum of the harmonic quantum oscillator (after a scaling 
and a shift). In particular, there is a spectral gap separating the origin from 
the rest of the spectrum. 

Finally, if one chooses a sequence {Pj^^} that converges to an appropriate 
limit sequence {Fj^^}, and if the initial n- velocity marginal in the limit N oo 
is given by 

/(")(«!, ...,v^-0) = fTivi, ^;„) + E E ^.-.^^^SH^i, ■■■,vn), (24) 

then its subsequent evolution is given by 

f^^\v,, ...,v^;t) = ft^{v,, ...,v^) + Y,Y. P^M^h^u • • • , «n)e"S^ (25) 
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and describes an exponentially fast approach to equilibrium in the infinite 
system. 

Coming now to the evolution equations for the /'^"'^(vi, v„; t), we here 
only state the final result, which is a special case of the more general one 
presented in the next subsection. In the limit ^ oo, the n-th marginal 
evolution equation for the diffusion master equation on M^^^^ = §^^^^, eq. 
()17|). becomes the (essentially linear) Fokker-Planck equation in M^", 

in particular, for n = 1 we recover eq. (0} for / = f^^\ with w = and Sq = e, 
and together with (jSJ and ((Tj). The momentum constraint © on the initial 
data is here immaterial. 

It thus appears that we have derived eq. (jU as a kinetic equation. Ap- 
pearances are, however, misleading. At this point Q does not yet have the 
status of a kinetic equation; notice that Kac's concept of propagation of chaos 
has not entered the picture! In fact, /*•"-' in p6|l may still in general be a 
convex linear ensemble superposition of extremal states, which are products of 
n one-particle functions representing the velocity density function of an actual 
member of the infinite ensemble. In other words, (j26p for n = 1, 2, . . . defines 
a "Fokker-Planck hierarchy" for a general statistical superpositions of initial 
conditions, analogous to the well-known Boltzmann hierarchy that arises in 
the kinetic theory of dilute gases |Spo91tl^IP94j . In this simple case, however, 
the n-th linear equation in the hierarchy (j26|) is decoupled from the equation 
for the n-l-l-th marginal. The upshot is that the first equation of the hierarchy 
is decoupled from /*^^) and therefore already a closed equation for f^^\ Since it 
is essentially linear (we say "essentially," for the parameter e is coupled with 
the initial data), a linear superposition of different solutions (corresponding 
to a statistical ensemble of initial data for / with same energy) is again a 
solution. Hence, if initially 

f^-\v,,...,v^;0) = {f^-{v,,...,v^)), (27) 
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where ( . ) is the Hewitt-Savage |HS55j ensemble decomposition measure on 
the space of initial velocity density functions of individual physical systems, 
where each fo{v) > satisfies (0) and ((Tj), then at later times 

/(")K,...,t;„;t) = (r>i,...,t;„;t)) (28) 

where f{v;t) solves eq. (jlj) with u = and Eq = s, and together with ^ 
and ((Zj). Note that the Hewitt-Savage measure is of course invariant under 
the evolution. This finally establishes the status of eq. (@]) (together with its 
constraints) as a kinetic equation valid for (almost) every individual member 
of the limiting ensemble. 

We remark that a product structure for f^^\vi, . . . , f„; 0) imposes inter- 
esting relations on the expansion coefficients, but we have no space to enter 
their discussion here. 

3.2 Gas in a periodic box 

A periodic box is physically unrealistic, but it provides a simple example of 
a situation in which momentum conservation has to be taken into account, 
too. In large parts the discussion of the previous subsection carries over to 
this situation. 

3.2.1 Finite 

The evolution of the ensemble of finite systems is now described by the 
diffusion equation on M^^""^, 

a,FW(V;t) = A„3.-4FW(V;t). (29) 

The spectrum of the Laplacian is now given by {A^i^^r_4 with ^^^3n~4 = 
j{j + 3N — 5), and the eigenspace for ^^^3N-i is spanned by the restrictions 

to c L^^^'^ of the harmonic polynomials which are homogeneous of 

degree j in L^^~'^; when j > this restriction to S^j^^^ C L^^~^ has to be 
non-constant. The computation of the eigenfunctions is again straightforward. 
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but the embedding of C L^^~^ causes a minor inconvenience because a 

rotation of all velocity variables is involved to get back to the physical velocity 
variables. We shall skip the details here, which will be presented elsewhere, 
and now turn directly to the derivation of the hierarchy of linear evolution 
equations for the marginal densities which obtains from eq. (j^HI)- 

It is convenient to express the Laplacian on the right-hand side in eq. 
fl29|) in terms of the projection operator from to the tangent space to the 
embedded manifold M^^~^. The relevant formula is discussed in the Appendix, 
eq. (fTTj). In order to apply eq. (f7T|) to eq. we introduce an orthogonal 

basis for the orthogonal complement in M.^^ of the tangent space to M^^~^ at 
V G M^^"*^. If li = 0, such a basis is simply provided by the set of vectors 
{V , El, E2, E3}, with E„ = (eo-, . . . , Co-), where the a = 1,2,3, are the 
standard unit vectors in M^. li u ^ 0, the E^r are mutually orthogonal but 
not orthogonal to V G M^^~^; however they are orthogonal to W = V — U. 
Indeed, 

1 ^ \ 

-J^E.^eA-V. (30) 



V-U={hN ^ 



(7=1 ^ 

The vectors {W , Ei, E2, E3} form the desired orthogonal basis for the or- 
thogonal complement in of the tangent space to M^^~^ at V G M^^~^. 
The magnitudes of the vectors E^ and W are \Ea-\ = VN and = ^/2Neo, 
respectively, where we recall that Eq = e — ^ \uf . Then, eq. becomes 

' 1 „A 



dV 



■ (31) 



(T=l 

In order to obtain an equation for the n-th marg inal of F(^)(V; t), which will 
be denoted by F^'^^^\vi, . . . , t), we integrate eq. (jHT|) over (f „+i, . . . ,vn) G 



p3N-3n 



Clearly, 

dv..,i ...dv, — .i-^ =J2 — . (32) 



Also, 



Q / 3 5ir(JV)\ 3 Q2p{n\N) 

dvn,. -dv.^-iY.E.^E.^UY^Y. (33) 

\cr=l / k=l i,j&n ■' 
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Finally, 



d ( (9FW\ 



dV \ dV J 



2- 9^-r^"'^-^^) + 



where tOj = Vi — u and the permutation symmetry of F^^^ was used. Putting 
everything together, if F^^^ satisfies the diffusion equation on M^^~^ it follows 
that the n-th marg inal satisfies 



3{N-n) ^ d 



2eoN 



3.2.2 The limit ^ 00 



The spectrum of — Ag3i v-4 in the limit iV — > cxd is given by 



lim I r = , (36) 



^0) 



which up to the replacement of e by £0 agrees with (j23|) . Thus, the whole 
spectrum is once again discrete and, in particular, there is a spectral gap 
separating the origin from the rest of the spectrum. 
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As a result, once again the velocity densities f^""^ approach equilibrium 
exponentially fast. The evolution equation for /*^") which obtains in the limit 
N ^ oo from (jHSj) is the (essentially linear) Fokker-Planck equation in M^", 



In particular, for n = 1 we recover eq. (@)) with / = f^^\ together with the 
constraints on the initial data (0), (©j and ((Tj). The last step to establish the 
status of (0)) (together with its constraints) as a kinetic equation involves once 
again the Hewitt-Savage decomposition, which implements Kac's concept of 
propagation of chaos for (|H7j). 

4 The Balescu-Prigogine master equation 

After having established that the linear Fokker-Planck equation (@]) together 
with the constraints on the initial data can be derived as a kinetic equation 
from eq. (j29|) . the diffusion equation on the energy-momentum foliation of M^^, 
we now study the more complicated diffusion process on the fohation M^^"'^ 
associated with the (nonlinear) Landau equation, eq. (P). At least at the 
formal level, the Kolmogorov equation for the diffusion process in question is 
given by the Balescu-Prigogine master equation |PB59j for the time evolution 
of F(^), which can be written once again as 




+ —{vi-u 



) 



(37) 



(38) 



where now 




(39) 



with 



v,w 



V — w 



(40) 
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and = lN\{k}; the density F^^^ has to satisfy the initial condition 

hm^j^o -^'•^'*(^; ^) = -^o^''(^)- ^ Fourier-transformed version of (j38p was con- 
structed in |PB57llPn59] . but it does not seem to have received much attention 
ever since. Balescu and Prigogine aheady pointed out that Landau's equa- 
tion (PJ can be extracted from (j38j) by contraction onto the first marginal of 
F^^\V;t). The formal argument runs as follows: Eq. (j38|) is equivalent to a 
hierarchy of evolution equations for all the marginals of F^^\V; t), which are 
denoted by t), with n = 1, A^. Of course, F^^l^) = F^^). 

The time evolution of F^^-^^^Vi, . . . , v^, t) is given by the corresponding con- 
traction of the Balescu-Prigogine master equation onto n variables, 

d,F^^\-^ = -^C^-^F^"\-^-^± [C.,,^„^^^^^^ (41) 

k=i 

Introducing the short-hands v for Vi and w for i>2, the evolution equation for 
the first marginal takes the form 

dtF^^\''\v;t)=d^-! \v-w\-^ Pi_^ - {d^-dJiF^^\'^\v,w-t)d^w. (42) 

Clearly, if we had F^'^\^\v,W]t) = F^'^^^\v;t)F'^^^^\w;t), then ^ would 
be a closed equation for F^^^^\ However, since (jHHj) with and (pUj) does 
not preserve a product structure of F^^^ for any finite A^, F^^'^^ cannot re- 
main a product of the first marginals even if that is the case initially. It was 
Kac's insight that if one lets ^ oo and assumes that limN^oo Fq^\v) = 
<S)kLi foi'^k), then one ought to be able to show that the product structure per- 
sists in time, viz. lim^too FW(y;t) = fi^k^t) for all t > 0, for which 
Kac coined the phrase "the propagation of molecular chaos." Kac proved 
propagation of chaos for his caricature of the Maxwellian gas, and Balescu 
and Prigogine surmised that propagation of chaos will hold also for their mas- 
ter equation when ^ cxd, in which limit equation (jl^ then becomes the 
Landau equation (^J). We remark that, for initial conditions that do not nec- 
essarily factorize, in the N —>■ oo limit eq. (pT|) leads (formally) to the Landau 
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hierarchy 

n „ 

9i/(")K,...,t;„;t) = -5^ / /:.„.„^j("+i)K,...,t;.+i;t)d^^„+i (43) 

k=i -^^^ 

A rigorous justification of propagation of chaos for the Balescu-Prigogine 
equation is an interesting open problem, on which we recently made some 
progress. Here we report on our results: 

1. the Balescu-Prigogine diffusion operator has an interesting geometric in- 
terpretation: it is a weighted average of the Laplacians associated with 
a certain family of sub-manifolds of M^^~^. These sub-manifolds are de- 
termined by the conservation laws for binary collisions between particles, 
as will be made clear below. 

2. for each finite the Balescu-Prigogine master equation is well-posed in 
£^ n £^ (M^^^'^) and displays exponential decay to equilibrium. 

3. as N —>■ oo, the first non-zero eigenvalue converges to zero. 

Our results show on the one hand that the finite N Balescu-Prigogine 
equation is more similar to the diffusion equation than meets the eye, yet on 
the other hand, the limit — > cxd is markedly different. In particular, our 
results suggest that the spectral gap vanishes. 

4.1 Geometric Aspects of the BP Master Equation 

Clearly eq. ()38|1 is a linear parabolic partial differential equation in 3A^ vari- 
ables, and it is easily verified that C^^\ like Aj^[3jv-4, annihilates the con- 
stant function m{V), the linear 3- vector polynomial p{V), and the quadratic 
polynomial e{V), ensuring the conservation of mass, energy and momentum. 
Hence, also in this case the evolution of the ensemble factors into independent 
evolutions on the invariant manifolds M^^~^. 

Now, for 1 < I < k < N consider the family of N{N —l)/2 two-dimensional 
manifolds 

®L = {V ■ Vk + vi = (Xki \vk -vi\^ = Pi Vi = j^l i^k,l^ (44) 
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where ctki, (3ki and ^'^^i are 3A^ — 2 arbitrary constants. Clearly, each point in 
M'^^ determines one such family of manifolds and each manifold corresponds 
to the conservation laws associated with the "collision" between particles k 
and /. In fact, the conservation oi \ vk — vi\ is equivalent to the conservation 



Vk + vi is also constant; thus, if V G ^ it follows that the manifolds B^^ 
determined by V are sub-manifolds of M^^""^ for all /c, /. 

Theorem 1. Let C^^w be the operator in eq. and F^^^ a probability density 
onMi^-^. Then 



- 
to 




\Vk\ ^ 





\vk-vi\-' AbjFW 



(45) 



where ^^1^ ihe Laplace- Beltrami operator on 

Proof. We first observe that the factor \vk — vi\^^ can be moved out of the dif- 
ferential operator in eq. (jlOl), since (^d^i — d^^) \vk — vi\^^ yields terms parallel 
to Vk — vi, which are annihilated by the projector P^^_^^. Hence, 

^v,,v, = -^\vk - vi\-^ ((9^, - 9^J-(P^,„„, ■ {d^, - ). (46) 

Next, we calculate ^^2^. At any point V G M.^^ the manifold B|; has the 
3A^ — 2 normal unit vectors M^j, E and Wi^a where 

/ \ / \ 



1 

71 











E = 



V^\vk - vi\ 





Vk - Vi 




-Vk + Vi 





{k entry) 



(47) 



(/ entry) 



V / V / 

and the e„, a = 1,2,3, are again the standard unit vectors in M'^; VFj^o- is the 
standard unit vector in M.^'^ with e„ as the i-th 3-block, i ^ k,l, and zeroes 
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everywhere else. Since all these vectors are mutually orthogonal, the projector 
on the tangent space to B^^ is 

P^2^ = hN - E ^ E -J^M^ ^ ^Wi,^ = (48) 





. -p^ 














k 


I 



k 
I 



where all the unmarked entries are zero and, of course, 

Vk-Vi Vk- Vi 

7 ® -: 



■"k-1>l 



\Vk-Vi\ \Vk-Vi\ 

It follows immediately from eq. (f7T| that 



(49) 

(50) 
□ 



which, together with eq. (jiUj) gives eq. (jl^ . 

Hence, the Balescu-Prigogine master equation can be written as 

= -T7^ E E K-«zr^ A«.F(^) (51) 



fceijv 



(fe) 



It is instructive to compare eq. (jHU), which leads formally to the Landau equa- 
tion, to the diffusion equation on M^^""^, eq. (UHl), which leads to the linear 
Fokker-Planck equation. In the "Landau" case the diffusion does not take place 
isotropically over the manifold of constant energy and momentum, but on the 
collection of the sub-manifolds B^^ determined by the two-particle conserva- 
tion laws. To be precise, the Balescu-Prigogine operator in eq. (j3T|) is a sort of 
weighted average of the Laplacians on all the B^^, where the fc/-Laplacian has 
weight ( "diffusivity" ) \vk — vi\^^ . These diffusivities are constant quantities 
on the corresponding manifolds, but they can take arbitrarily small values in 
certain regions of as increases. Loosely speaking, the ellipticity of the 
BP equation degenerates as N ^ oo. This is bound to affect the rate of decay 
to equilibrium of the solutions in the same limit. 
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4.2 Well-posedness of the BP master equation and de- 
cay to equilibrium 

We now identify the evolution of F^^\V;t) with the motion of a point ipt 
in the Hilbert space £^ (M^^""^) . We define the Sobolev-type space as the 
closure of C°°(M^^~^) w.r.t. the norm || . ||^^, given by 

11^11?/= Q^''nV^,V^) + ll^lli(M--), (52) 

where Q^^^ is the manifestly symmetric positive semi-definite quadratic form 
associated with the operator C^'^^ defined on (M^^~^) : 

gW(^/;,0)= / ^jJC(^)(j)dT (53) 

for {tp, 4>) in (M^^-^) The form closure of Q^^\ip, 0) in ^ x i^, which is 
also denoted by Q^^^{ip, 0), defines a unique self-adjoint, positive semi-definite 
operator with dense domain ^ C ^, the Friedrichs extension of C^'^\ also 
denoted by C^^\ We recall that the Sobolev-type space 9j coincides with the 
domain of the square root of the Friedrichs extension of C^^\ i.e. the operator 
; see |RS80j for general background material. Thus C^^'^ is a densely 
defined, unbounded operator on the Hilbert space £^ (M^^~'^) . 

It is easily seen that the kernel space of Ker(£*^^)) = DIq is one-dimensional, 
consisting of the constant functions. Hence, £^ (M^^~^) decomposes as DIq © 
where is the orthogonal complement of DTq in £^(M^^~^), i.e. the 
subspace of £^(M^^~'^) on which C^'^^ is strictly positive. Also, 9j decomposes 
as © we remark that Sj^ C ii^'"*" can be equivalently defined as the 
closure of G C°°(M3^-^) : J^sN-i'^dr = 0} w.r.t. the norm || . , where 



^,^='g(^)(^,V^). (54) 

Since C^^^ is self-adjoint and strictly positive on i^+fl^, it follows immediately 
that the operator C^^'> is the generator of the contraction semi-group e^*^'^' 
on |RS75j . This implies that C^^^ is also the generator of a strongly 



21 



continuous semi-group on all of S? (M3^-4) , also denoted by e"*^^"^' . Precisely, 
e"*^'^'' acts isometrically on the invariant subspace D^o and strictly contracting 
on its orthogonal complement hence e~*^^^' is also positivity preserving. 
Thus, if G 2?{ml^-^), then 

V;f)^='e-*^'"Vr (55) 

solves uniquely for Cauchy data limfjo V'*-^-' = V'o^^ and the initial value 
problem for the BP master equation is well-posed. The evolution of an initial 
density 4""^ G H ^'(M^^-^) (i.e., 4"^^ > and /m3^^-4 V^^^dr = 1), 
actually takes place in the intersection of the positive cone of 2? (M^^~^) with 
the affine Hilbert space 

2l = V'(^)+£2'+, (56) 

where t/^^^ = jM^^-^p^ is the Hilbert space vector of the (analog of the) 
micro-canonical equilibrium ensemble. 

The spectrum of can be studied with the standard techniques de- 
veloped for weak solutions of linear inhomogeneous PDE in divergence form 
|GT98j . extended |Heb96j to operators on compact manifolds without bound- 
ary (here M^^~^). Some care must be taken due to the fact that the ellipticity 
is not uniform, since the coefficients are unbounded above; on the other hand, 
the ellipticity condition is satisfied uniformly from below. First of all, since the 
bilinear form Q^^^ is clearly positive and bounded on ij^, by the Lax-Milgram 
Theorem it follows that for A < the operator c'^^'^ = C^^^ + AX determines 
a bijective mapping from S)^ to {Sj^)* (the dual of Sj^). Next, we introduce a 
compact embedding S from 9)^ to (i?^)* such that {Su){v) = J^3n-4UV dr for 
all V G S)^. In order to define the operator S, we first observe that points on 
M3^-4 satisfy 

\vk - vi\^ = \vkf + \vif - 2vk ■ vi < 2{\vkf + \vif) < me (57) 
so that 

Q(^\^|J,^|J)>^=J- V C^(^,^) (58) 
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where 



= 11 {d^, - d^,)ij ■ P^^^^^ ■ K -9.>dr. (59) 

Since M^^~^ is a compact manifold, it is easy to see that there is some constant 
Cat such that 



(60) 



where Vip is the covariant derivative of iJj on M^^~^ and jV"?/^!^ = g"^^ diip djip , 
g"^^ being the metric tensor and djtp the derivative with respect to the j-th 
coordinate. Combined with eq. (j3Hj) this gives 

II^IlL > , ||^||L,2 (61) 

"^'"^ 2y2iVF(iV-l)"^"^+ ^ ^ 

where 2J^^ is the closure of {ip G ^""(M^^^^^ : J^3N-4t/jdT = 0} w.r.t. the 
norm 



2 dci 



■3JV-4 



iVipf dr. (62) 



Equation (jHT|) implies that i^"*" is continuously embedded in 211^^, which in 
turn is compactly embedded in by the Sobolev embedding theorem and 
the Rellich-Kondrasov theorem, both of which hold on a compact manifold 
|Heb96j . Finally, ii^'"*" is continuously embedded in (S)^)* (via the Riesz Rep- 
resentation Theorem), and the compact embedding S of i^"'' into (^^)* is 
obtained as the composition 5)+ —>■ W^^"^ £?'^ —>■ (^^)*. 

Finally, one obtains the standard Fredholm alternative |GT98j . by re- 
writing the equation C^^^tp = a as 

^ + (A-Ao)<)^^ = <V (63) 

where ^^^^ (Aq < 0) is the (continuous) inverse of jC^;^\ Then, since S is 
compact, —(A — Xo)&^^^^S is also compact from S^^ to Sj^, and the standard 
Riesz-Schauder theory of compact operators in a Hilbert space leads to the 
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conclusion that C^^^ has a purely discrete spectrum and that each eigenvalue 
has a finite-dimensional eigenspace. Since the spectrum of C^^^ is discrete, 
we have a spectral gap between Aq^^ = and the smallest non-zero eigenvalue 
A^^^ > of C''^\ and we conclude that the equilibrium ensemble is approached 
exponentially in time 



(N) 



£2 



(JV) 



{N) 



■ (64) 



4.3 The limit iV ^ oo 



The question whether the spectral gap remains finite in the limit — > cx), 
which was recently answered affirmatively in the context of the Kac model and 
other models related to the Boltzmann equation |CCL02j , and also by us in the 
previous section for the diffusion master equation, has presumably a negative 
answer for the BP master equation, at least for the Coulomb case studied here. 
Indeed, we now show that the smallest positive eigenvalue with permutation 
symmetric eigenfunction vanishes in the limit, i.e. lim^v^oo -^i^'' = 0. 
Consider 



(TV) def 



inf (V',>C(^V) 



(65) 



where (using permutation symmetry) 



N 
~2 



[ {d., - J V' • I • {d., -d.,)tp dr (66) 



Here, S+ = G Sj-, 



3JV-4N = 1}, and S)'^ is the permutation 



symmetric subspace of Sj^. An upper bound on a[^^ is obtained by selecting 

We use 



a suitable trial function ip G 



with g{v) 
C = 



N 

1=1 



-c 



(67) 



v1/2. In order to satisfy the condition if) G one has to choose 



A 




?>N -I 



M37V-4 



(68) 
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Then, the quantity ( 'th, C^^'ib ) can be calculated exactly. In the 

"standard" case ^ u = 0,5 = 1, a tedious exercise leads to the estimate 

Ai^) < ^ ^ (iV > 1) (69) 

Thus, A[^'' ^ as ^ oo. 

This is not enough to conclude rigorously that the BP equation has a van- 
ishing spectral gap a,s N ^ oo (because we have not determined the asymp- 
totics for all the eigenvalues). However, in view of our previous remarks on 
the fact that the ellipticity of the equation degenerates as — > oo, it seems 
reasonable to conjecture that this may well be the case. In turn, a vanishing 
spectral gap for the BP equation would lend support to the conjecture |Vil98aj 
that the Landau equation itself possesses solutions that decay to equilibrium 
slower than exponentially. 

We end on the remark that this result can be easily generalized to any 
BP master equation corresponding to a Landau kinetic equation with "soft" 
potentials (for 7 < 7/3 in eq. Q). Conversely, it is not difficult to prove that 
for BP master equations with "hard" potentials (7 > 7/3), the spectral gap 
remains finite in the limit N —>■ 00. 
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5 Appendix 

We here recall some basic facts about the Laplacian on a Riemannian manifold. 

Let M be a n-dimensional Riemannian manifold with metric Qij, and let / be a 

^In general, the evolution of Pj;^'' on M^^"** can always be obtained from the evolution 
on Mq^~^ via the simple transformation V U + eoV. 
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function on M, / G C°^(M). The familiar expression for the Laplace-Beltrami 
operator acting on / in local coordinates is 



where g = det{gij). Now, let us take M to be embedded in the Euclidean space 
M"^, m > n, and let Qij be the metric induced by the standard metric in M*". Let 
the function / be defined and differentiable over some suitable sub-domain of 
M™ containing M. We associate to any point a; G M a set of orthonormal basis 
vectors ei, . . . , e„, e„+i, . . . , such that Ci, . . . , e„ span the tangent plane 
T-cM, whereas e„+i,...,em span the orthogonal complement T^M.-^ in M*". 
The expression for A/ is coordinate- independent, and we are free to choose 
local coordinates (xi, . . . such that the coordinate curves are tangent to 
the Ci, . . . , e„. In such coordinates gij = 6ij, g = 1 and d^^f = (V/)-,-, i.e the 
component along ej of the gradient V/ of / in M™. Then eq. ()70|) becomes 



where Pm is the orthogonal projection from to T^M.. Note that eq. (jTlj) 
is coordinate-independent in R'", and preserves permutation symmetry with 
respect to the Cartesian coordinates of M"*. This is important when dealing 
with probability densities for A^-particle systems, which are constrained by the 
dynamical conservation laws to evolve on certain lower-dimensional manifolds 
and at the same time must be permutation-symmetric. 
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